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Abstract 

Using confocal microscopy we investigate the effect of attraction on the 
packing of polydisperse emulsions under gravity. We find that the distri- 
butions of neighbors, coordination number, and local packing fraction as a 
function of attraction are captured by recently proposed geometrical mod- 
eling and statistical mechanics approaches to granular matter. This extends 
the range of applicability of these tools to polydisperse, attractive jammed 
packings. Furthermore, the dependence of packing density and average coor- 
dination number on the strength of attraction provides the first experimental 
test of a phase diagram of athermal jammed particles. The success of these 
theoretical frameworks in describing a new class of systems gives support to 
the much-debated statistical physics of jammed matter. 

The random packing of particles is a problem of great practical and theoretical 
importance with applications to systems as diverse as granular materials, emul- 
sions and colloids. How athermal particles pack is affected by many parameters, 
including the packing protocol, the shape and roughness of the particles, as well 
as their polydispersity. Nevertheless, experiments and numerical simulations have 
suggested the existence of reproducible disordered packing structures, which can 
be achieved by sedimenting and shaking monodisperse ll25l [l5l and polydisperse 
ifTTTl spheres, ellipsoids ETI and tetrahedra [16]. For frictional particles, there 
emerges a reversible packing density curve from random loose packing (RLP) to 
random close packing (RCP) as a function of the applied tapping protocol f22\ . 
The dependence of density of these jammed states on the friction coefficient has 
led to an equation of state and a prediction of a jamming phase diagram to de- 
scribe the behavior of these systems in numerical simulations |[26ll . The repro- 
ducibility of all these jammed states for a given set of experimental conditions 
suggests that a statistical mechanics framework may apply to these inherently out 
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of equilibrium systems [ 13 ]. Other studies have argued that the history dependence 
and heterogeneities in the packing preclude a general theory of jammed granular 
matter HU ED. 

In this paper, we introduce a novel route to jamming by packing athermal, fric- 
tionless spheres under gravity in the presence of a short-range attractive force. By 
examining the packing geometry in 3D we find that the strength of the attractive in- 
teraction controls the packing density and the microstructure of jammed states. We 
thus provide the first experimental verification of another branch of the phase di- 
agram of jammed matter originally predicted for monodisperse repulsive particles 
|[26ll . namely that of polydisperse attractive particles. Moreover, our experimen- 
tal distributions of local packing properties are in very good agreement with the 
predictions of recently proposed models applied to repulsive particle aggregates: 
the 'granocentric' model [11] and the 'i-gamma distribution' of the local packing 
fraction @. This study therefore extends the limits of applicability of geometri- 
cal modeling and statistical mechanics approaches to capture the effects of both 
attraction and polydispersity. These results open the way to use tools originally 
developed to investigate granular materials to a new class of systems. 

Jamming attractive athermal particles under gravity should not be confused 
with the widely studied gel transition in thermal systems, where classical thermo- 
dynamics applies [7]|20l|24]]. In thermal systems attractive forces lead to colloidal 
aggregates with fractal structures at arbitrarily low densities, while here the absence 
of thermal fluctuations seems to localize the effect of attraction to the first shell of 
neighbors where the forces due to gravity and attraction balance to form a mechan- 
ically stable pack. Despite the fundamental differences in the mechanism of their 
creation and the resulting structures, this study reveals one important similarity. 
We observe a reentrant behavior of packing density with attraction, analogous to 
that observed in thermal glassy systems ll24l l27l . thus questioning its postulated 
dynamical origin. 

To create a jammed packing we cream athermal oil-in-water emulsions under 
gravity, as described in reference [9 1. In addition, we introduce an attractive poten- 
tial by varying the concentration of sodium dodecyl sulfate (SDS) micelles, which 
act as a depletant [6]. The depletion energy has an entropic origin and its de- 
pendence on the surfactant concentration is derived in reference [3] and validated 
for our emulsion system in ||6l |U [T4j [T2J. We therefore calculate the short-range 
attractive force F^ at an equilibrium distance of around 10 nm and using a 
critical micellar concentration (CMC) of 13 mM [10|. The range of attraction is 
limited by the saturation point of SDS (i.e. 50 mM) at a maximum average force 
of (Fa) = 33 pN. This average depletion force is defined for each packing as the 
force between two droplets of average radius, (r) = 3.5 /urn. Since the particle size 
distribution has a 25% standard deviation in radius, the depletion force between 
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Figure 1 : Confocal 2D slices and their reconstructed images for repulsive (a, b) and 
attractive (c, d) 3D emulsions. Ticks mark the original and reconstructed spheres. 
Navigation maps (white lines) are superimposed with the packing (b, d). Voids 
shown in (c) appear to occupy cells as particles above and below the confocal 
plane share void space in the third dimension. 

particle pairs picked at random, calculated using the standard AO theory |]3j , has 
a distribution with ±20% standard deviation. This calculation is justified because 
the the micelles are approximately 500 times smaller than the size of the small- 
est emulsion droplet and the micelle volume fraction is always below 2% lfl4l[T2ll . 
The average attractive forces (F^) studied here range from 10 to 100 times stronger 
than the weight of a single particle. This corresponds to small deformations in the 
range of 1 — 10 A, such that the spherical approximation holds. 

The protocol to prepare the jammed packing involves creaming a very dilute 
emulsion at a density (j) 5% to avoid aggregation of the particles. Refractive index 
matching the droplets and the aqueous phase allows us to image the dynamics of 
the packing process using a fast scanning confocal microscope (Leica TCS SP5 
II). The voxel size is 130 nm in the horizontal xy plane and 300 nm in the vertical 
z direction. We image a box size of 65 x 65 x 100/um, several particles away from 
the boundaries of a \mL container. Once the packing of « 1500 droplets is formed, 
confocal images in Figs. QJ, c reveal a visual difference between the attractive and 
repulsive jammed packings: repulsive packings display homogeneous structures 
while high enough attraction gives rise to voids within the packing. We analyse the 
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structure in terms of the particle positions and radii with subvoxel accuracy using 
a Fourier transform algorithm [0. We then use the geometrical overlaps in the 
reconstructed spheres (see Figs. [TJ), d) to identify particles that are in contact. This 
allows us to measure the number of contacts a particle has, i.e. its coordination 
number, z. Error in estimating particle positions and radii translates to an error of 
±0.3 in the (z) estimation. To characterize the local neighborhood of each particle 
we tessellate space using the navigation map El [TTJ shown in Figs. \Vp, d. This 
mapping defines two local quantities: the number of neighbors, n, a particle shares 
an interface with and the local packing fraction, (j)^ oc , which is the ratio between 
the particle volume, Vp, and the volume of the corresponding cell, V. The global 
density, (j), is determined by the total volume of the particles divided by the system 
volume. 

In Fig. 12 we present the probability density distributions of the local parameters 
Z, n, and (j)j oc as a function of attraction. Even though the packings exhibit more 
void space as the attraction is increased, the neighbor number distribution P(n) 
shown in Fig. [2^, is independent of attraction. Just as in the repulsive packings, 
the (n) is 14 ±0.3 with a standard deviation of 30%. This makes sense because the 
navigation map partitions the extra void space between the neighbors, but the num- 
ber of common interfaces remains the same. The inset shows that the particle size 
dependence on the number of neighbors also remains the same for all attractions, 
rendering the number of neighbors an insensitive measure of the packing structure. 
The coordination number z reports on the number of force-bearing particles in the 
mechanically stable packing, which is greatly affected by the attractive forces in 
addition to gravity. The measured distributions P(z) in Fig. (2J) peak at decreasing 
values from 7 to 4 and become narrower as a function of attraction [19j]. The iso- 
static condition predicts a value of (z) = 6 for repulsive frictionless spheres in 3D 
CQ, yet the inset in Fig. [2j) shows that (z) decays as a function of (Fa) to below 
this value. With increasing attraction the probability of such local configurations 
also increases, as shown in the P(z), thus shifting (z). Analogously, the presence of 
stable loose voids shifts the distributions of local packing fraction ^((|>i oc ) shown 
in Fig. [2h towards lower packing densities with wider distributions as a function of 
increasing attraction. The decrease in global density (j) is also apparent in the top 
panel of Fig. |2ji. 

On the other hand, the initial increase in (z) above isostaticity and (j) above RCP 
at the onset of attraction, shown in the graphs to be significantly larger than the er- 
ror bars on the measurements, remains puzzling. While jammed random states 
above RCP have been proposed numerically in athermal systems ll28l . this is the 
first experimental system to explore such compact configurations that avoid crys- 
tallization. Whether these denser states become accessible due to polydispersity 
remains an open question. This turnover implies that isostatic packings with a den- 
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Figure 2: (Color online) Probability density distributions P(n) in (a), P(z) in (b), 
and ^(^loc) in (c) are shown as a function of (F^), while the insets in (a, b) present 
the dependency of the average values. The legend in (a) defines the colors and 
symbols for all data sets, which are used throughout the manuscript. Repulsive 
emulsion packings at RCP are included for comparison. All distributions are fitted 
by the granocentric model, whose parameters are shown in the inset in (c). The 
global trends in (F^) and (z) as a function of (j) are shown in (d). 
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Figure 3: (Color online) (a) Evolution of P(V) and its cumulative distribution for 
different levels of attraction. Notice that V is rescaled by the volume of the average 
particle for each distribution, (b) Data collapse is observed when V is rescaled by 
the shifted average volume of the cell (V). The inset shows the compactivity % 
decrease with increasing density (j). 



sity at RCP can also be achieved with a nonzero attractive potential, as shown by 
the magenta points and those at 19 pN in both panels in Fig. |2ji. Interestingly, a 
similar turnover in versus (j) is seen in attractive colloidal gels undergoing struc- 
tural arrest and interpreted as the signature of a reentrant glass transition Il24ll27l . 

To understand the observed trends in the probability distributions in Fig. [2] 
quantitatively, we interpret the data using the granocentric local model introduced 
in ifTTl . The model uses as inputs the experimental values of (n), (z) (excluding 
particles with z < 3) and global density (j), and outputs the distance 8g between two 
neighboring sphere surfaces, as well as the distributions P(n), P(z), and P((j)j oc ), 
shown in Fig. |2k,b,c. In the inset in Fig. |2fc, the ratio of contacts to neighbors, 
given by p = ((z) — 3)/((n) — 3) is shown to decrease as a function of the attractive 
force due to the decrease in (z) , while the predicted distance 8g naturally increases 
as the packings become looser and agrees very well with the average experimental 
value (8exp)- Further confidence in the model is shown in Figs. [2^, b, c by the 
agreement between the measured distributions and the theoretical predictions in 
which all parameters are fixed by experimentally determined values. It is surpris- 
ing that a local stochastic model captures all the probability distributions, since it 
assumes that there are no local correlations between neighboring particles nor long 
range correlations beyond the first shell of neighbors. This result is in contrast 
to the structure of attractive colloidal gels, where the fractal dimension indicates 
correlations between thermal particles that persist throughout the system 11711201. 

Given the success of the local granocentric description of the packing structure, 
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we next consider the fluctuations in the volume V of each cell in the navigation map 
in terms of a statistical mechanics framework [13 1. Independently picking volumes 
V from a uniform distribution results in an exponential probability distribution of 
volumes P(V) in the thermodynamic limit. This Boltzmann-type distribution max- 
imizes the entropy given the constraints. Assuming that each cell V is made of k 
elementary cells (5]|3, this distribution becomes a shifted &-gamma distribution 
with a shape parameter k and a scale parameter (V) — V imn : 



Indeed, a wide range of monodisperse packings with different packing protocols 
and global densities have been successfully fit by this type of distribution ll5ll4"l[T7l 
[181 . Here we investigate the effect of polydispersity and subsequently attraction 
on P(V). Fig. [3^ shows that polydispersity of 25% alone significantly broadens the 
distribution compared to the previously measured monodisperse case (star sym- 
bols) [5 ] and shifts the global packing fraction to a higher value. Rescaling each 
V by the volume of the particle itself leads to a broader distribution /^loc) in the 
polydisperse case, shown in Fig. [2};, indicative of an additional source of random- 
ness to the size distribution. Fitting the mono and polydisperse cases with k rj 14 
and k « 2, respectively, captures the observed changes in shape and width. This 
means that the polydisperse packing approaches the Boltzmann distribution. As a 
function of attraction, Fig.[3k shows an increase in both {V} and the standard devi- 
ation a of the distributions, as expected from the looser packing structures. Since 
Eq. (Q]) gives ((V) — V imn )/a = y/k, this dependence leads to a constant value 
of k = 2.0 ±0.4 for all levels of attraction. This is shown in Fig. |3j) by the data 
collapse that results from rescaling the volumes as (V — V m j ri )/((y) — V imn ). 

Within the thermodynamic framework, a simple counting argument of volumes 
leads to a definition of the entropy S and consequently an analogue of inverse tem- 
perature = dS/dVi i, which is known as the compactivity in granular statistical 
mechanics. The resulting expression for % = ((V) — V rrun )/^ is shown to be lin- 
early decreasing with density in the inset in Fig. [3j>. Recently, a phase diagram 
for jammed matter have been proposed in terms of the dependence of % and (z) 
on density for monodisperse frictional hard spheres ll26l . In comparison, our data 
shown in Fig. [2]i lies below the predicted RLP line at infinite % (dash-point line), 
which we have extended beyond monodisperse RCP for comparison. We achieve 
packings denser than RLP over the same range of (z) , consistent with the fact that 
we explore states with decreasing values of %. This also shows that the line we 
observe is not simply the RLP line shifted due to polydispersity. Surprisingly, we 
extend the limits of the phase diagram to hyperstatic packings up to (z) = 7.5 and 
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densities up to 74%, which is above the RCP density (68% in the poly disperse 
case). The microscopic mechanism for jamming with attraction and the origin of 
these states provide interesting questions for further study. 
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